Entanglement renormalization for quantum fields 
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It is shown how to construct renormalization group flows of quantum field theories in real space, as 
opposed to the usual Wilsonian approach in momentum space. This is achieved by generalizing the multi- 
scale entanglement renormalization ansatz to continuum theories. The variational class of wavefunctions 
arising from this RG flow are translation invariant and exhibit an entropy-area law. We illustrate the con- 
struction for a free non-relativistic boson model, and argue that the full power of the construction should 
emerge in the case of interacting theories. 
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Classical statistical mechanics, quantum many-body sys- 
tems, and relativistic quantum field theories all involve an 
extremely large number of degrees of freedom living at dif- 
ferent length scales. The interactions between these de- 
grees of freedom are the source of notorious difficulties in 
their study. However, much insight has been gained from 
the renormalization group (RG), which has proven to be 
the natural tool to deal with the different length scales in 
such systems 01. In its original development, the RG acts 
as a fixed operation at the level of the classical partition 
function (or related quantities such as the effective action). 
This operation is typically formulated in momentum space 
and can only be implemented exactly for free theories. Per- 
turbative expansions in a small parameter are required for 
interacting theories. In addition, this formulation of the RG 
is only applicable to quantum systems using the quantum- 
to-classical mapping, which is known to fail in some cases 

a. 

One exception is Wilson's numerical renormalization 
group (To, which can be interpreted as an implementa- 
tion of the RG directly at the level of the quantum wave 
function. Together with White's more powerful density 
matrix renormalization group (DMRG) Q], these meth- 
ods are now understood as a variational optimization over 
the class of matrix product states (MPS) @]. Based on 
the observation of an entropy/area law JD] in locally inter- 
acting quantum lattices, quantum- information-theoretical 
considerations have resulted in the development of other 
sophisticated variational classes for quantum lattice sys- 
tems. These are generally known as tensor network states 
and can be associated with RG schemes, allowing the clas- 
sification of gapped phases of matter f6[]. Unlike Wilson's 
fixed RG scheme, these are variable RG schemes that can 
be variationally optimized. They are formulated in real- 
space and deal equally well with free and interacting sys- 
tems. One specific scheme, called entanglement renormal- 
ization 0], can also be applied to study critical phases and 
can be used to compute, e.g., scaling exponents Hfl. The 
corresponding variational class, the multiscale entangle- 
ment renormalization ansatz (MERA) is set apart by 
its unique properties, including, the ability to support al- 



gebraically decaying correlations and logarithmic correc- 
tions to the entropy/area law in (1 + 1) dimensions. This 
class has been successfully applied to study the physics of 
a wide variety of strongly interacting systems in low di- 
mensions, including the study of real-time evolution, and 
fermionic and anyonic systems which are inaccessible by 
Monte Carlo methods J^-Ell. 

Most of these developments have been restricted to the 
lattice setting. While they do allow the study of contin- 
uous quantum systems via discretisation, it is often de- 
sirable to work directly in the continuum. Examples in 
condensed matter physics include strongly interactin g ul - 
tracold atomic gases lfl3ll and impurity problems [Ql 14], 
whereas the fermion doubling problem IU7II clearly moti- 
vates a continuum treatment of relativistic theories. Re- 
cently it was understood how to define MPS and its higher- 
dimensional analogue for continuous systems II 1 511 . This 
class has already been successfully applied to study both 
non-relativistic and relativistic quantum field theories [[l6ll . 
It is the objective of this Letter to define a generalisation of 
MERA directly in the continuum. The area law for entan- 
glement entropy — with logarithmic violations for critical 
theories in (1 + 1) dimensions — has also been observed in 
the continuum in the context of black hole physics iflHl and 
conformal field theories lfl9ll . and thus validates the poten- 
tial usefulness of our approach. 

The MERA class. — The MERA construction, introduced 
in |@|) may be described either via an active renormal- 
ization process applied to a strongly correlated quantum 
state or, dually, as the result of a special quantum cir- 
cuit applied to a simple fiducial state. This reverse de- 
scription is not in violation with the irreversibility of the 
RG, since it only applies to the ground state, not to the 
whole Hilbert space and the Hamiltonian acting upon it. 
At stage of the MERA construction, the degrees of free- 
dom (distinguishable quantum spins in this case) are ar- 
ranged in a regular lattice and initialised in some conve- 
nient initial state, e.g., the "all Os" state |0). At stage 1 
the degrees of freedom are subjected to a local interaction 
Ui for some constant time, resulting in a correlated quan- 
tum state Ui |0) of the lattice. The precise details of XJ\ — 
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while playing a key role in applications of the variational 
principle — are not required to establish the general proper- 
ties of a MERA. At the next stage the lattice is subjected 
to a scale transformation and the lattice spacing is dou- 
bled. At this point, in order to renormalize the lattice (i.e., 
restore the lattice spacing), new degrees of freedom are in- 
troduced: one adds a quantum spin initialised in the state 
|0) between each pair of the old lattice sites. We write this 
renormalization step as 1Z. The resulting state IZUi |0) 
is then again subjected to a local interaction U 2 for some 
constant time followed by the renormalization step 7Z, pro- 
ducing the state IZU^RUi jO). This process is then iter- 
ated as many times as desired/required. A physical in- 
terpretation of the MERA construction is straightforward. 
Every layer k begins with a dilation of the state living in 
the previous layer to a lattice doubled in size, by applying 
1Z. The operation of Uf. then adds short-range (e.g. over 
two sites) fluctuations/entanglement on top of this state, 
corresponding to fluctuations over 2 m ~ fc sites in the final 
state if m layers are applied. However, the resulting state 
I^mera) = U m 1UJ m - 1 'Tl- ■ - IZUi |0) generally breaks 
translation invariance unless the Uj are carefully fine tuned. 

The passage to the continuum. — For simplicity we spe- 
cialise to the case of a single bosonic species in one di- 
mension. (The generalisation to higher dimensions and 

to fermions or Bose-Fermi mixtures is entirely straightfor- 
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only degrees of freedom which are nontrivially affected are 
those with momenta k < A. 

The third ingredient of the MERA construction is the 
renormalization step where the scale is changed and new 
uncorrected degrees of freedom are introduced. Here the 
continuum analogue is not entirely clear, but we argue 
that the following replacement naturally achieves the same 
objective: we simply effect a small change of scale via 
W = exp(— iJL), where the generator L is given by 

x , V O) dx. (2) 
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This has the same physical outcome as the original MERA 
renormalization step because initially uncorrected degrees 
of freedom at lengthscales below e are now introduced at 
the new lengthscale. Thus our proposal for the contin- 
uous MERA (cMERA) is as follows: evolve some ini- 
tial state \Q) according to K for an infinitesimal time <5, 
which correlates real-space degrees of freedom at scales of 
0(e). Then introduce new degrees of freedom from the 
higher momenta/shorter lengthscales by dilating the state 
via evolution according to L for a time 5. The last layer 
s = s e creates fluctuations at lengthscale e. No fluctua- 
tions at a shorter range exist in the final state. The fluctu- 
ations created by layer s live at lengthscale e exp(s e — s) 
in the final state. If this process is to correctly produce 
the long-range entanglement in a state with correlation 



ward.) We write ^x) and ^ (x), x £l, for the bare field length ^ ^ first , ayej . shouM bg defined at g = with 



annihilation and creation operators which obey the canoni 
cal commutation relations [i/>(x), ^{y)\ = S(x — y). The 
following constructions can also be described in terms of 
any set of operators that define the theory, such as the her- 
mitian field operator (f> and its conjugate momentum ir for 
relativistic boson theories. 

The generalisation of the Oth stage of the MERA con- 
struction is clear: one should choose for the initial state a 
factorized reference state The transition to a continu- 
ous space x £ R enables the introduction of a continuous 
scale parameter s that labels the layers of the MERA con- 
struction. In every layer, new fluctuations are created by 
the action of a unitary evolution U(s) = exp(—i5K(s)) 
with time step 5 and local interaction K(s) given by 



(1) 



K(s) = / k(x, s) dx. 



where k(x,s) is a local combination of the field opera- 
tors ift(x), their adjoints, and their derivatives. Since the 
continuum lacks a shortest distance, we need to introduce 
a characteristic lengthscale e below which K(s) does not 
create entanglement. This is possible in a variety of ways: 
one strategy is to build k(x) from smoothed field opera- 
tors ip e (%) = J Xe(x — y)ip(y) dy, where Xe( x ) is some 
smooth envelope function which is nonzero outside a re- 
gion of width e around x = 0. Another strategy is to sim- 
ply impose a cutoff on K when it is expressed in terms of 



momentum variables. The cutoff A 



ensures that the 



s e — s 5 = 0(log(£/e)). By taking the limit 5 — >■ we 
obtain our final expression: 

where T is the time-ordering operation. We refer to the 
unitary operation preparing a cMERA as U(s e ,s^). Note 
that a UV cutoff A = e _1 and an IR cutoff are explic- 
itly built into the cMERA definition. For critical systems 
or relativistic theories with an infinite range of fluctuations 
it is necessary for s e — — > oo. 

The biggest difference between the cMERA and MERA 
definitions arises from the flexibility we have in imposing 
the UV cutoff. In the MERA case the UV cutoff is dic- 
tated by the lattice spacing. This, in turn, essentially forces 
the form of the subsequent scaling transformation (i.e., an 
integral number of spins must be introduced in the scaling 
step). The freedom, arising from the continuum, to choose 
a smooth UV cutoff allows the scaling transformation to be 
applied continuously and for the resulting state to be easily 
chosen to be translation invariant. 

The set of all cMERA forms a variational class: the 
variational parameters are the coefficients of interactions 
involved in K; these coefficients may depend on both x 
and the integration parameter s. When the coefficients do 
not depend on position x a generic cMERA is manifestly 
translation invariant (this is established by noting that ap- 
plication of the unitary e lSL to a translation-invariant state 
results in a translation invariant state). This is in contradis- 
tinction to the generic situation with lattice MERA. 
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Comparison with Wilsonian renormalization. — The 
cMERA has been constructed using the quantum circuit in- 
terpretation, but can now be interpreted as an active renor- 
malization process and compared to Wilson's momentum- 
shell renormalization group (RG) [1]. The latter works by 
iteratively integrating out all the modes living in a small 
momentum shell A — dA < \k\ < A. A rescaling step 
brings modes at A — dA back to A. In the end A can be 
sent to infinity, but we need to start with a finite A in order 
to define the process. The renormalization process defined 
by the cMERA proposal is a real-space implementation of 
Wilson's momentum shell RG in a hamiltonian framework: 
rather than integrating out high-frequency modes around 
the cutoff A from the partition function — which is a fixed 
operation — the operator K first disentangles these modes 
from the wavefunction in such a way that they can be iso- 
metrically projected onto a reference vacuum |f2). Then 
a scale transformation is performed to send the disentan- 
gled modes beyond the cutoff (where they no longer in- 
teract via K) and new entangled modes are brought to A. 
These modes are then disentangled in the subsequent step. 
This immediately clarifies the need for a finite cutoff A in 
our construction (which can also be sent to infinity at the 
end of the process). 

So what are the differences between the renormalization 
process defined by the cMERA and Wilson's momentum- 
shell RG? Firstly, whereas Wilson's RG is a fixed opera- 
tion, the cMERA renormalization process is governed by 
K{s) which can be variationally optimized. Secondly, 
while K{s) can be formulated in momentum space {e.g. 
to implement the cutoff), its defining property is real-space 
locality, which has proven to be a correct assumption for 
both free and interacting theories in countless examples 
with MERA and related variational classes and is a result 
of the locality of physical interactions. 

The cMERA RG flow. — New to our formalism is that 
we can define the RG flow for operators in a hamilto- 
nian framework: suppose we want to evaluate the expec- 
tation value (O) = (*|0|^), where |^} is a cMERA 
of the form Eq. d3j, and O is some local operator, e.g. 
O = ip(0) or O = ip' (O)ifj(x). To do this we define 
O(s) = U(s e , S yOU{s e , s), where 



^l = -i[K(s)+L,0(s)]. 



(4) 



This 'heisenberg-like' equation of motion is obtained by 
differentiating the lower limit of the evolution operator 
U(s e , s). Then (O) may be found by integrating this equa- 
tion from s = s e down to s = s^, with the initial condition 
0(s e ) = O, and evaluating (fi|0(s 5 )|n). 

Physically, we think of the bare or physical operator O 
as being defined at the UV cutoff lengthscale x ~ e. As 
"time" s progresses we think of O(s) as being brought 
from lengthscale x ~ e to lengthscale x ~ ee s '~ s . Thus, 
O(s) is obtained from the the bare operator O by renor- 
malizing up to scale s (i.e. all degrees of freedom between 



momentum scales e s ~ Sc A and A, where s < s e , have been 
integrated out/disentangled). 

For critical theories, K(s) is expected to become s- 
independent away from s e . In accordance with the results 
from JD], we can then assume the existence of operators O 
that satisfy — i[K + L, O] = AO. These are scaling op- 
erators with scaling dimension A. If O is a local scaling 
operator, it is necessarily centered around x = 0. A local 
scaling operator 0{x) with scaling dimension A at position 
x satisfies 

-i[K + L,0] =xdO(x)/dx + XO(x). (5) 

The existence of scaling operators makes it easy to prove 
that cMERA support algebraically decaying correlations 
and are thus well suited to study critical theories. In ad- 
dition, we can illustrate that they are able to produce an 
area law for the entanglement entropy. 

An entropy/area law for cMERA. — We now provide a 
heuristic argument that a generic cMERA obeys an en- 
tropy/area law by appealing to results [l20ll concerning the 
dynamics of quantum spin systems: it is known that the 
entanglement entropy SA(t) = — tv(pA(t)log(pA(t)), 
where pA(t) = tiA"{^ ltH |<A)) ('/'ol e ltH ) is the reduced 
density operator for a region A, under the real-time dynam- 
ics generated by any strongly interacting system grows as 



dS A (t)/dt < c\dA\ 



(6) 



where \dA\ denotes the length or area of the boundary of 
A and should be measured in terms of the cutoff of H, for 
some constant c which depends only on the local interac- 
tions and the geometry of A. It is natural to conjecture that 
this result holds in the continuum setting. Now, subject 
to this assumption, we can deduce the proposed area law: 
we track the entropy growth of the time-dependent region 
A(s) = Ae s ~ s % i.e. A(s) is A scaled down by a factor 
e s-s e throughout the cMERA preparation. We bound the 
entropy Sa by integrating Eq. © (compare with iI^Hi ): 



S A <c (LAe 

Js,-log(LA) 

fclog(LA), 
M^LA)- 1 (l 



s — s € \d—l 



ds 



(LA)- 



d= 1 
d> 1, 



(7) 



where the area \dA\ < (LA) d_1 ; the entanglement of 
A(s) with the remainder of the system cannot receive fur- 
ther contributions when |A(s)| < {LA) d < 1. The ap- 
pearance of the logarithm of the UV cutoff is familiar from 
standard QFT calculations OGl]. Note that the cMERA 
might not describe logarithmic violations of the boundary 
law in d > 1, similar to the MERA case. 

Representation of ground states via cMERA. — Let us 
now construct a cMERA representation of the ground state 
of a simple non-relativistic bosonic model 



dx. 
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which is well-defined if 2v < fi. For 2v = [i, the el- 
ementary excitation becomes massless and the model be- 
comes critical. Using a general strategy discussed in [I22II . 
we can describe ground states of free theories with a Gaus- 
sian cMERA where K(s) is the quadratic operator 



K(s) 



fg(±s) fi(k)fi(-k)-$(-k)$(k) 



dk, 



2s + A 2 /A 4 e- 2s ) + 2/i/A 2 ] , 
and A = (/i 2 - 41/ 2 ) 1 / 2 



where ip(k) is the Fourier transform of the field opera- 
tors. We set g(k/A,s) = x{s)T(\k\/A) with T(k) a 
fixed cutoff function with cutoff 1. The variational param- 
eters are thus the function values xi s )- 1^) is fixed by 
= 0. The analytic calculations are facilitated using 
a sharp cutoff such as T(k) = 9(1 — \k\), where 6(x) is 
the Heaviside function. Although this cutoff function pro- 
duces a nonlocal operator K, it is straightforward to see 
that similar results are obtained using a smooth cutoff such 
as r(/c) = exp( — k 2 ) which does yield a local K. The 
exact ground state of H can be accurately reproduced if 
A 2 3> O(n) by choosing 

X (s) = 2^/A 2 )/ [(e 

where we have set s e 
represents the mass gap of the system. For 2u < fx or 
thus A > 0, the 'disentangling strength' x( s ) peaks at 
s = — l/21og(A 2 /A) and decays to zero for s —> — oo. 
The integration of the RG flow can be stopped at <C 
-l/21og(A 2 /A) - C(log(^/A 2 )). In the critical limit 
(2u — > fi), x( s ) reaches a non-zero horizontal asymptote 
lim s _ ) ._ 00 x( s ) = v l ' I 1 = 1/2 and we need to integrate all 
the way down to = — oo. According to Eq. (01, the scal- 
ing operators then correspond to <f)(x) ~ [ip{x) + ip^(x)) 
and ir(x) ~ (ip(x) — ijj\x)). The low-energy behavior is 
scale-invariant and can be described by the massless Klein- 
Gordon equation rf22ll . 

Conclusions. — In this Letter a generalisation, cMERA, 
of the MERA variational class to the continuum setting 
has been introduced. We have argued that cMERA can 
be translation invariant and generic ally exhibit an en- 
tropy/area law. We have also supplied an analytic argument 
that the ground states of a general class of local quadratic 
models admits a cMERA description. Much remains to be 
done: we expect, by analogy with MERA, that cMERA 
will be a useful variational class for strongly interacting 
quantum fields, and will allow the description of a vari- 
ety of interesting physical phenomena from topological ef- 
fects to confinement, and symmetry breaking. Looking 
further afield, the cMERA constitutes a realization of the 
holographic principle. It is tempting to speculate, build- 
ing on [23] and 12411 , that cMERA are a natural candidate 
to establish a link between entanglement renormalization 
and the best known realization of the holographic princi- 
ple, namely the AdS/CFT correspondence. 
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